Coulomb blockade and Kondo effect in the electronic structure of Hubbard molecules 
connected to metallic leads: a finite-temperature exact-diagonalization study 
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The electronic structure of small Hubbard molecules coupled between two non-interacting semi- 
infinite leads is studied in the low bias- volt age limit. To calculate the finite-temperature Green's 
function of the system, each lead is simulated by a small cluster, so that the problem is reduced to 
that of a finite-size system comprising the molecule and clusters on both sides. The Hamiltonian 
parameters of the lead clusters are chosen such that their embedding potentials coincide with those 
of the semi-infinite leads on Matsubara frequencies. Exact diagonalization is used to evaluate the 
effect of Coulomb correlations on the electronic properties of the molecule at finite temperature. 
Depending on key Hamiltonian parameters, such as Coulomb repulsion, one-electron hopping within 
the molecule, and hybridization between molecule and leads, the molecular self-energy is shown to 
exhibit Fermi-liquid behavior or deviations associated with finite low-energy scattering rates. The 
method is shown to be sufficiently accurate to describe the formation of Kondo resonances inside the 
correlation-induced pseudogaps, except in the limit of extremely low temperatures. These results 
demonstrate how the system can be tuned between the Coulomb blockade and Kondo regimes. 

PACS numbers: 73.23.Hk, 73.21. La, 72.15.Qm,+a, 73.20.At 



I. INTRODUCTION 

Finite-size electron systems linked to non-interacting 
electron reservoirs have been a topic of intense theo- 
retical and experimental study because of their rele- 
vance to quantum dot systems and single-molecule de- 
vices. Depending on the importance of correlation ef- 
fects induced by the electron-electron Coulomb repul- 
sion, different types of theoretical approaches are em- 
ployed. For weakly correlated systems, ballistic elec- 
tron transport is studied within the one-electron approx- 
imation such as density- functional theory (DFT)i~— On 
the other hand, strongly correlated systems are mod- 
eled by tight-binding Hamiltonians with Hubbard- or 
Anderson-type interaction terms and various many-body 
techniques are applied £~— Two noticeable effects be- 
yond the one-electron approximation are the Kondo ef- 
fect and Coulomb blockade, both of which are observed 
in quantum dot systems^ - — More recently, the Kondo 
effect was also observed in adsorbed molecules by scan- 
ning tunneling spectroscopy and high-resolution photoe- 



mission spectroscopy. 
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The Kondo effect in nano-size 



systems was studied theoretically by using the numeri- 
cal renormalization group (NRG) technique for a variety 
of cases, such as multi-dot or multi- level systems and 
dots coupled to superconducting leads ^r— within DFT 
combined with the one-crossing approximation— - — and 
the continuous-time quantum Monte Carlo technique. 27 
Coulomb blockade effects seen in electron transport 
through a finite-size interacting system are investigated 
by using rate-equation techniques and non-equilibrium 
Green's function theory^ - — 

The aim of the present work is to introduce a new 
scheme for the investigation of quantum dots that is ap- 



plicable in the full range between Kondo physics and 
Coulomb blockade, except in the limit of extremely low 
Kondo temperatures. To illustrate this approach, we 
focus on small interacting molecules coupled to non- 
interacting semi-infinite electron reservoirs. The many- 
body properties of these systems are evaluated by using 
exact diagonalization (ED) at finite temperatures^"— 
In order to apply ED, the semi-infinite leads are simu- 
lated by finite-size clusters. For a given chemical poten- 
tial, the tight-binding Hamiltonian parameters of these 
clusters are chosen such that the difference between the 
surface-site Green's function of a semi-infinite lead and 
the corresponding cluster lead is minimized along the 
Matsubara axis. The finite-temperature Green's func- 
tion of the total system consisting of molecule and lead 
clusters is then evaluated exactly within ED. Since the 
effective lead-cluster Hamiltonian is extremely sparse, at 
typical temperatures of interest only a limited number 
of excited states needs to be evaluated. Here we con- 
sider Hubbard chain and ring molecules attached to two 
metallic leads. In these systems lead clusters consisting 
of only five bath levels can accurately mimic the true 
embedding potentials down to T > ^m/500, where tM 
denotes the hopping parameter representing the metallic 
leads. Thus, the formation of Kondo resonances within 
the correlation-induced pseudogaps can be investigated. 
In the case of single leads, the problem reduces to the 
single-impurity Anderson model. The cluster size can 
then be significantly increased so that much lower tem- 
peratures can be reached. An additional advantage of 
our method is that the zero-bias- voltage limit of the in- 
teracting system can be studied for arbitrary values of 
the Coulomb repulsion and molecule-lead hybridization, 
without further approximations. 
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The discretization of the semi-infinite leads and the 
application of ED to the resulting finite-size system is 
analogous to the use of ED as impurity solver within 
the context of dynamical mean-field theory (DMFT) and 
its cluster extensions^"— where the Weiss mean-field is 
represented by a finite number of non-interacting lev- 
els. Multi-orbital as well as multi-site correlations have 
been studied for a variety of materials, 43 including vari- 
ous surfaces and heterostructures4£ The main difference 
is that in the present case the leads are assumed to be 
non-interacting, so that the self-consistent iterative pro- 
cedure is absent. 

An important feature of the ED approach is that it 
provides complete dynamical information, in particular, 
transfer of spectral weight between low and high excita- 
tion energies, formation of Hubbard bands, and opening 
of correlation- induced pseudogaps. A quantity of cen- 
tral interest therefore is the molecular self-energy which 
exhibits strong variations as a function of Hamiltonian 
parameters and temperature. In particular, Fermi-liquid 
behavior in the Kondo regime and correlation-induced fi- 
nite scattering rates in the Coulomb blockade regime can 
clearly be identified. 

The outline of this paper is as follows. In Section HT1 we 
describe our theoretical model for the molecule-lead sys- 
tem and discuss several details of the calculation of the 
molecular self-energy and interacting Green's function. 
In Section IIIII we present the numerical results and the 
discussion, with special emphasis on the Coulomb block- 
ade and Kondo effect. Section HVl contains the summary. 
In the Appendix we consider a single adatom on a semi- 
infinite lead and examine the temperature range in which 
a finite-size cluster can be used to simulate a semi-infinite 
lead. Since this case is equivalent to the single-impurity 
Anderson model, the results can directly be compared 
with those of other schemes that are applicable at arbi- 
trarily low temperatures. 



II. THEORY 

A. Formalism 

We consider a molecule consisting of N atomic sites 
and linked to two semi-infinite metal leads, as shown 
schematically in Fig.QJa). The isolated molecule is mod- 
eled by a single-site Hubbard Hamiltonian characterized 
by the onsite energy e a , the nearest-neighbor hopping 
interaction t, and the onsite Coulomb repulsion U, 



k = E« 
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(hj)cr 



cLci 



h it hii, (1) 



where c\ a (c^) creates (destroys) an electron with spin a 
at site i (1 < i < A/"), hi a = c^c^, and the summation in 
the second term is taken over pairs of nearest-neighbors. 
Hereafter, we adopt the notation where the matrix (oper- 
ator) corresponding to a quantity A is denoted by A while 




semi-infinite cubic lattice M-site cluster 



FIG. 1: (a) Tight-binding model for an N-site molecule at- 
tached to two non-interacting semi-infinite leads, (b) For the 
evaluation of the Green's function of the molecule, the two 
semi-infinite leads are replaced by a cluster consisting of M 
levels, as shown here schematically for the left lead. 



its matrix elements are defined as Aij . For simplicity, we 
limit the discussion here to one level per molecular site 
and purely onsite Coulomb interactions. Equivalent ly, it 
would also be feasible to investigate multi-orbital inter- 
actions, including inter-orbital Coulomb and Hund's rule 
coupling, for instance, in transition metal ions attached 
to semi-infinite leads. Throughout this paper, the hop- 
ping integral within the molecule is taken as unit of the 
energy scale, i.e., t = 1. 

The left (right) lead is represented by non-interacting 
electrons on a semi-infinite simple cubic lattice with 
nearest-neighbor hopping interaction tM and the onsite 
energy level is chosen as zero of the energy scale: 



h 



L(R) 



J2 tM 



(2) 



where z, j < 1 (i, j > iV), so that the energy bands of 
both leads exhibit a finite density of states (DOS) in the 
energy range [-6£ M ,6£ M ]- 

The molecule is linked to the left (right) lead via the 
hopping integral between site 1 (N) of the molecule and 
site (7V + 1) of the left (right) lead. The mixing term of 
the Hamiltonian describing the molecule-lead hybridiza- 
tion is expressed as 



hm.ix — 



t R c 
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AT+lcr 



CNc 



h.c, 



(3) 



where, for simplicity, the hopping integrals on both sides 
are assumed to coincide: tL = £r. The Hamiltonian of 
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the total system consisting of the molecule and the two 
leads is given by 



H = he + h L + h R + h Tl 



(4) 



We investigate the electronic structure of the molecule 
described by Eq. (j4j) for a wide range of Hamiltonian pa- 
rameters. As we consider the low bias- volt age limit, both 
leads have the same chemical potential which is denoted 
as fi. The Green's function G of the molecule can be 
written as 



Gij(iuj n ) = iuj n + fj,- he - H(iu) n ) 

- s L (iuj n + /i) - s R (iuj n + /i) 



(5) 



where 1 < z, j < TV, u n = (2n + 1)ttT (n > 0) are 
Matsubara frequencies at temperature T, h^ denotes 
the first two terms of Eq. (pQ), and E(zco> n ) is the self- 
energy matrix accounting for electron correlation effects 
within the molecule. As we also allow for interatomic 
Coulomb correlations, the self-energy matrix has off- 
diagonal components with respect to site index. This 
approach differs from recent ones22r— in which the self- 
energy of each site is assumed to be local and determined 
in a self-consistent manner similarly to the layer DMFT 
approach^ In the present work, we consider only param- 
agnetic solutions and omit the spin index a hereafter. In 
Eq. ([5]), s 1 "^ denotes the embedding potential describ- 
ing the one-electron hybridization effects due to the left 
(right) lead on the molecule j 49 i 5Q These embedding po- 
tentials give rise to broadening and shifting of the molec- 
ular levels and therefore play the role of contact self- 
energies^ 

For the present geometry, only the {11} element of s L 
is non- vanishing: 



*n(*) = *i 



00 



= t\ 9oo(z) ■ 



(6) 



Similarly, the only non-vanishing element of s R is 
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The surface Green's functions appearing in these expres- 
sions are given by 



6tA 



de 



(8) 



where pi(e) denotes the local density of states per spin at 
surface site i within the left (i = 0) or right (i — N + 1) 
lead. 

In order to make use of ED for the evaluation of the 
interacting Green's function of the molecule, Eq. (j5j), we 
follow a procedure that has proved to be very useful in 
analogous DMFT calculations. The surface Green's func- 
tions gu(z) (i = 0,7V + 1) representing the continuous 



spectra of the leads are approximated by those of finite 
clusters consisting of M levels, as depicted in Fig. djb) 
for the left lead. The {00} element of the cluster Green's 
function is given by 



g c o l o(z) 



e 



M-l 



k=i 



(9) 



where e k (0 < k < M - 1) and v k (1 < * < M - 1) 
are the energy levels and intra-cluster hybridizations, re- 
spectively. An analogous expression holds for the surface 
Green's function of the right lead, ^tv+i,at+i(^). (Note 
that the cluster levels e& do not refer to actual lattice 
sites within the leads. Instead, they represent auxiliary 
quantities to simulate the spectral distributions of the 
leads.) The discretization of goo(^) is not suitable on 
the real energy axis since ^oo(^) has a continuous energy 
spectrum while #oo( z ) possesses only a finite number of 
poles. Thus, lmgQ l (z) — > or — oo in the limit z —> 0, 
whereas goo(^) for metallic leads remains finite, so that, 
in the very low-energy region, the cluster Green's func- 
tion deviates strongly from the actual lead Green's func- 
tion. These discrepancies are absent if the calculation is 
restricted to finite temperatures, goo(^) can then accu- 
rately be fitted by g^z) at Matsubara frequencies, since 
both functions vary smoothly along the imaginary energy 
axis. As shown in Ref. 43, in finite-T ED/DMFT calcu- 
lations for typical multi-orbital materials, two or three 
bath levels per orbital are adequate to achieve adequate 
fits for temperatures in the range T w VK/50, . . . , W/200, 
where W is the bandwidth. For the present case, this im- 
plies T w 0.025, . . . , 0.10. To reach lower temperatures 
therefore requires accordingly larger lead clusters. 

As in standard ED/DMFT calculations, for a given 
chemical potential /i, the discretization of goo(^) can be 
achieved by determining e& and v k in Eq. (|9]) via mini- 
mization of the quantity^ 
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] W n | goo(iw n + AO - 9oo( iuj n + aO 



(10) 



where the weight function W n is chosen as l/u n in or- 
der to provide greater accuracy at low uo n . (The large 
frequency behavior is less relevant in this fit since both 
Green's functions approach l/{iuj n ) at large uo n .) Other 
choices, such as W n = 1 or W n = usually give very 

similar results, even though the auxiliary cluster param- 
eters €k and Vk may vary slightly. With decreasing tem- 
perature T the lowest Matsubara frequency approaches 
the real energy axis, so that the fitting becomes less ac- 
curate. As will be demonstrated in the next section and 
Appendix, the true lead Green's function can be simu- 
lated by that of a relatively small 5-level cluster with 
sufficient accuracy as long as the temperature is approx- 
imately in the range T > £m/800. Moreover, the lower 
boundary of this temperature range can be reduced by 
increasing the cluster size. As a consequence, it is feasi- 
ble to describe the Kondo effect on the spectral density, 
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if the associated Kondo temperature is comparable with 
this temperature range. 

It should be noted here that the Matsubara tempera- 
ture used in the fitting of the lead surface Green's func- 
tion may be viewed as a fictitious temperature Tm that 
does not need to coincide with the physical temperature 
T. Instead, its choice is mainly determined by the num- 
ber of cluster levels used to simulate the semi-infinite 
leads. Evidently, a larger value of M permits fitting at 
lower values of Tm- This point will be addressed further 
in the Appendix where an extremely small value of Tm 
is chosen for the evaluation of the self-energy of a single 
adatom over a wide range of real temperatures. In most 
applications discussed below the Matsubara temperature 
is taken to be the physical temperature. 

Let us denote the non-interacting Green's function of 
the molecule linked to the two semi-infinite leads by G° 
and that linked to the two clusters by G°' d , where the 
term 'non-interacting' signifies U = in the molecule. 
The interacting counterparts are G and G cl , respectively. 
When the tight-binding parameters of the lead clus- 
ters are optimized as described above, one can presume 
that G°(iu n ) ~ G 0,cl (iuj n ) at all Matsubara frequen- 
cies. As a result, when the Coulomb interaction in the 
molecule is switched on, G(iuo n ) should nearly coincide 
with G cl (iuj n ). We may therefore employ G cl as a reason- 
able representation of the true interacting Green's func- 
tion, G, of the molecule attached to the two semi-infinite 
leads. Schematically, the procedure outlined above pro- 
ceeds via the following steps: 
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Below we do not distinguish between the molecular 
Green's functions G and G cl . We emphasize, however, 
that even if G agrees well with G cl at Matsubara points, 
at real energies G cl has a discrete level spectrum while 
that of G is continuous. 

The Hamiltonian of the interacting molecule linked to 
the two M- level clusters is highly sparse. To compute 
G cl (iu n ), we therefore make use of the Arnoldi algo- 
rithm which is ideally suited to evaluate the lowest eigen 
states relevant at temperature T. G cl {iuo n ) is then de- 
rived via the Lanczos procedure for a finite number of 
excited states. In the present work, system sizes up to 
n s = N + 2M = 15 have been investigated. More details 
concerning the numerical procedure are provided in Refs. 
|42| and l43l. (As mentioned earlier, the approach outlined 
above does not involve any self-consistency procedure as 
in DMFT since the leads are uncorrelated. Their fixed 
electronic structure merely governs the boundary con- 
ditions of the interacting molecule. The auxiliary lead 
quantities e& and Vk characterizing the leads are there- 
fore determined before carrying out the exact diagonal- 
ization.) 
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FIG. 2: (Color online) Comparison of the {00} element of the 
Green's function of a semi-infinite lead (solid circles) and that 
of a 5-level cluster (crosses), (a) ji = 0, (b) fi = 1, and (c) 
ji = 2, for t M — \ and T = 0.02. The real part of goo vanishes 
for \i — 0. 



B. Semi-infinite vs. cluster leads 

As discussed above, the calculation of the electronic 
structure of an interacting molecule between semi-infinite 
leads is made feasible by simulating the {00} ({iV+1, N+ 
1}) surface element of the Green's function of the left 
(right) lead in terms of a cluster Green's function, as 
indicated in Eq. ((9]). To demonstrate the accuracy of this 
fitting procedure, we compare in Fig. [2] both quantities 
as a function of uo n for three values of fi for a semi-infinite 
lead with t M = 1 (band width W = 12). 

The lead clusters consist of five levels. Thus, for \i ^ 0, 
there are in total nine independent fit parameters: e& 
(0 < k < 4) and v k (1 < k < 4). At half-filling (/x = 0), 
this number is reduced to four because of symmetry rea- 
sons: eo = and the other four levels are symmetrically 
distributed with respect to e = 0. The fitting then be- 
comes slightly less accurate than away from half-filling. 
For the present choice of T = 0.02, goo(iu n ) is seen to 
agree very well with gQ l (iuj n ) in the whole u n range. At 
small Matsubara frequencies weak cusps appear in the 
cluster Green's function as a result of its singular behav- 
ior along the real energy axis. 

Figure [2] suggests that the embedding potentials of the 
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ear and ring molecules, where the true leads are replaced 
by clusters consisting of five levels. 



C. Spectral information 

To demonstrate how the molecular electronic structure 
undergoes a transition between the Kondo and Coulomb 
blockade regimes, we consider in the next section the par- 
tially integrated quasiparticle density of states which is 
defined as 



PiOi) = -G«(r = 0/2), 



(12) 



where Gu (r) is the diagonal component of the imaginary- 
time Green's function at site i (/3 = 1/T): 
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(13) 



FIG. 3: (Color online) Non-interacting Green's function of a 
linear molecule with N = 3 between two semi- infinite leads 
(solid circles) and that between two 5-level clusters (crosses), 
(a) Diagonal {11} element and (b) off-diagonal {13} element, 
for t M = l,t L = 0.5, n 1.0, and T 0.02. 



semi-infinite leads in Eq. (j5]) can be approximated at 
Matsubara points by those of small clusters. To illus- 
trate this point further, we compare in Fig. [3] the resul- 
tant non-interacting Green's function of a linear molecule 
(N = 3) between two semi-infinite leads, G^-(io; n ), with 
the one of the same molecule between two 5-level clusters, 
Gij Cl (iuj n ). The Hamiltonian parameters correspond to 
those in Fig. [2fb) and the contact integrals are chosen 
as ti, = 0.5. It is seen that both the diagonal and off- 
diagonal elements of G^{iuo n ) are in excellent agreement 

with the corresponding cluster elements G^ cl (iuj n ). 

A crucial question determining the usefulness of the 
approach outlined above concerns the range of temper- 
atures in which accurate results can be obtained for a 
given cluster size. To explore this point, we present in 
the Appendix a careful study of the electronic structure 
of a single correlated adatom on a semi-infinite lead. In 
this case, a large range of cluster sizes can be employed 
in order to systematically investigate the behavior of the 
self-energy at very low temperatures. Cluster sizes up 
to M = 11 were used for T > £ M /1600. The results 
for 5-site clusters are found to agree quantitatively with 
those of larger clusters for T > £m/500, and qualita- 
tively for T > £m/800. These results, together with the 
ones shown in Figs. [2] and [3j demonstrate the usefulness 
of our strategy of evaluating the Green's function of the 
Hubbard molecule by simulating the semi-infinite metal- 
lic leads in terms of finite clusters. In the following, we 
present results for the electronic structure of various lin- 



and the interacting quasiparticle DOS is defined by 

1 



-lmGu(uj + iS), 



(14) 



with a positive infinitesimal 5. pi(p) may therefore be 
expressed as 

Pi(») = -J2e~ iu ^ /2 G u (iLo n ) 

n 

/oo 
doj F(oj) pi(u), (15) 
-OO 

where F(u) = 0.5/ cosh(/^/2) = (-Tdf/duj) 1 ^ 2 is a dis- 
tribution of halfwidth w = 41n(2 + y/S)T = 5.268T cen- 
tered about uj = (/ is the Fermi function). Thus, at 
a given value of the chemical potential /i, pi{n) repre- 
sents the quasiparticle DOS of site i partially integrated 
within a few T around \i. The advantage of studying 
this quantity, as compared to the actual interacting DOS 
pi(cj), is that it can be evaluated without extrapolating 
the Green's function from Matsubara frequencies toward 
the real energy axis. 

We also note that the weight function F(lj) in Pi(p) = 
—Gu(f3/2) is closely related to the one appearing in the 
conductance of single atoms attached to leads. In this 
case the weight function is given by 51 F c (uo) = —df/doo 
with halfwidth w c = 21n(3 + y/2)T = 3.525T. Thus, 
Pi(p) samples the molecular density of states near p in 
a window ~ 1.5 times larger than in the case of the con- 
ductance. 

We point out here that one can, of course, also evalu- 
ate the molecular Green's function Gf, close to the real- 
energy axis. However, these spectra consist of many 
sharp peaks related to the finite number of levels of the 
lead clusters. Moreover, the level energies e& and intra- 
cluster coupling terms Vk depend not only on the cluster 
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size, but also on the choice of the weight function W n 
in Eq. ([TP]) . Evidently, one would have to use very large 
lead clusters so that these discrete spectra evolve into a 
meaningful representation of the continuous spectra of 
the actual semi-infinite leads. Because of the exponen- 
tially growing Hilbert space the ED approach would then 
no longer be practical. Thus, the purpose of introducing 
the auxiliary fit parameters e& and Vk is to simulate the 
finite-temperature lead Green's functions and embedding 
potentials at Matsubara points. At not too low T, these 
functions converge very well with cluster size and are re- 
markably stable against variations in and and for 
different choices of W n . A2> The continuous spectra of the 
Green's function elements Gij(u) describing the molecule 
attached to semi-infinite leads may then be derived via 
analytic continuation of Gij(iuj n ) or, preferably, of the 
molecular self-energy T>ij(ioj n ) to real energies. In the 
latter case, continuation of the known one-electron prop- 
erties of the molecule and of the leads is avoided. The 
task of analytically continuing Gij{iuo n ) or T>ij(ioj n ) is 
entirely analogous to the one in quantum Monte Carlo 
simulations, where the maximum entropy scheme is of- 
ten used to generate real energy spectra. In DMFT stud- 
ies the discrete Green's functions Gfj(uj) can, however, be 
very useful for the identification of a Mott transition since 
the excitation gap opens at the same critical Coulomb in- 
teraction as in the continuous spectrum of Gij(u). 



(a) 

' t =1, t =0.5 

M ' L 

U=0 






J 









2 


(b) 
















3 


J 







-4 -2 2 4 

8 



FIG. 4: Non-interacting local DOS averaged over all sites for 
(a) linear and (b) ring molecules (N = 4) between two semi- 
infinite leads. The molecular configurations are depicted in 
the insets. £m = 1, £l = 0.5, and U = 0. Imaginary part of 
energy, 7 = 0.01 



III. RESULTS AND DISCUSSION 

A. Non-interacting molecules 

To provide an impression of the electronic structure 
of the molecule in the absence of Coulomb interactions, 
we show first in Fig. \3\pav( e ), the non- interacting local 
DOS averaged over all sites, for linear and ring molecules 
with N = 4 linked between two semi-infinite leads, where 
tM = 1 and ti, = 0.5. In the following, the onsite en- 
ergy e a is specified as — E//2, so that the system becomes 
electron-hole symmetric when it is half-filled. For lin- 
ear molecules, the DOS consists of TV resonant peaks 
corresponding to the energy levels e m of the isolated 
molecule, which are distributed symmetrically with re- 
spect to e = 0. The lowest energy state has even parity 
with respect to the center of gravity of the molecule, and 
the parity alternates in the order of ascending energy 
levels. The energy width of the m-th resonance, which is 
determined by its coupling to the imaginary part of the 
lead embedding potentials, is given by 

T m = 7Tt 2 L [po(em)|^m(l)| 2 + pN+1 (e m ) (N) | 2 ] , 

(16) 

where i/j m (i) denotes the amplitude at site i of the elec- 
tron wave function of the m-th level of the isolated 
molecule with energy e m . It is understood that the lowest 
level corresponds to index m = 1. For the ring molecule 
in Fig. H(b), the DOS exhibits only three peaks, since 



the central one at e = is doubly degenerate. The wave 
functions of the isolated ring molecule in the site basis 
are, 
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among which ^2 and ^3 are degenerate at e = 0. Since 
ip2 has no amplitude on sites 1 and 4 which are coupled 
to the leads, ^2 remains a truly localized interface state, 
even when the coupling to the leads is introduced, so that 
it makes a S- function contribution to p av (e). To avoid 
this singularity, an artificial imaginary energy 7 = 0.01 
is used in Fig. |4] for this level. The widths of the other 
levels correspond to the physical broadening. It should 
be noted that, once U is switched on, ^2 and ips are 
mixed, so that both states become delocalized. 



B. Coulomb blockade 

It should be emphasized that, in contrast to bulk sys- 
tems, where the ratio U/t (or U/W) provides a measure 
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FIG. 5: (Color online) (a) Partially integrated quasiparticle 
DOS, pi(aO — —Gu(/3/2), of a linear molecule with N = 3 
as a function of chemical potential \i for £l = 0.3, 0.5 and 
1.0 (U = 2, t M = 1, and T = 0.02). Small numbers near the 
minima of these curves indicate chemical potentials giving 
integer occupations, (b) Electron occupation for the same 
parameter set as in panel (a). 



of the strength of electron correlations, in the present 
case an important parameter characterizing the elec- 
tronic structure of the interacting molecule in the vicinity 
of molecular level e m is the ratio U/T m , where T m is the 
level width defined in Eq. (PS). For U/T m <C 1, the 
molecule is expected to be in the ballistic regime where 
correlation effects are dominated by the molecule-lead 
hybridization. On the other hand, for U/T m ^> 1, the 
molecule is in the Coulomb blockade regime where the 
onsite Coulomb repulsion hinders the addition of the sec- 
ond electron to the molecule when the first one occupies 
a resonant level. 

Figure[5](a) shows pi(3)(/i) = — ^11(33) (P/2) for a chain 
molecule (TV = 3) as a function of /i for three val- 
ues of the contact integral The other parameters, 
t m = 1, U = 2, and T — 0.02, are common to all curves. 
The integrated quasiparticle DOS of site 2, /?2(a0 = 
— G22OS/2), ( n °t shown) is similar to — Gn(/3/2), except 
that the peak structure arising from the second energy 
level around \i = is absent. Since this molecular or- 
bital has odd parity with respect to the center of the 
molecule it has no weight on site 2 irrespective of the 



magnitude of U. For tL = 1, intra- molecular correla- 
tion effects are dominated by single-particle hybridiza- 
tion with the leads. Thus, pi(/i) = — Gn(/3/2) exhibits 
three broad peaks as a function of /i, which originate 
from the energy levels of the non- interacting molecule. 
With decreasing II (increasing U/T m , see Eq. (fTE])). the 
DOS peaks start exhibiting minima at their centers. For 
tL = 0.5 and 0.3, all three peaks split into pairs of peaks 
separated by a pseudogap induced by Coulomb blockade. 
Note that, upon decreasing the double peaks on both 
sides of the minima become sharper, while the energy 
separation between them depends only weakly on tL- 

To analyze this trend, it is useful to expand the 
Coulomb repulsion in Eq. (pQ) in terms of the orbital basis 
of the isolated molecule. Specifically, the density-density 
interaction components are given by 

m,n 

where h^ m(T denotes the occupation of the m-th orbital 
with spin a. For the linear molecule with N = 3, one has 
U n = U33 = 3/7/8, U22 = U/2, U 12 = U 23 = U/4, and 
Uis — 3/7/8. The effect of other non-diagonal elements 
not included in Eq. ([18]) is small if U/t is not large. In 
the limit of small tL, the unrestricted self-consistent field 
(USCF) approximation may then be used to estimate the 
mean-field values of the m-th molecular level with spin 
a: 

&rri(j — 6 m -|- ^ ^ Ujyui (^"0 n ,— <j) • (19) 
n 

Thus, the first molecular level yields peaks at —^t — U/2 
and -y/2t-U/8, with a gap 3/7/8. The second level has 
peaks at — Z7/4 and Z7/4, with a gap U/2, while those of 
the third level are located at y/2t + U/S and y/2t + U/2, 
with a gap 3/7/8. The energy positions of the DOS peaks 
in Fig. [5ja) are seen to be in fair agreement with these 
estimates. 

Figure E^b) shows the total electron occupation of the 
molecule, summed over spin and site components, as a 
function of fi. For t L = 1, the occupation varies smoothly 
from zero to six. In contrast, for tL = 0.3 it exhibits 
distinct plateaus at each integer occupation, where those 
corresponding to odd integers are caused by the Coulomb 
blockade effect and their energy positions coincide with 
those of the pseudogaps shown in panel (a). 

To illustrate the effect of Coulomb correlations in the 
vicinity of the pseudogap, we plot in Fig. [6] the imag- 
inary part of the diagonal element of the self-energy, 
ImEn (icu n ), for a chain molecule with N = 3. The 
three curves correspond to tL = 1, 0.5, and 0.3. The 
other parameters are the same as in Fig. [5l Panel (a) 
shows the self-energy at half- filling (/i = 0). For tL = 1, 
the system is Fermi- liquid- like since ImEn(za; n ) tends to 
zero linearly as uj n 0. This behavior is in accord with 
the shape of the corresponding partially integrated DOS, 
Pi{v) — —G\\{fi/2), shown in Fig. [5fa), which exhibits 
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FIG. 6: (Color online) Imaginary part of correlation-induced 
self-energy, — ImEn (ioo n ), for a linear molecule with N = 3 
at (a) fi — (half- filling) and (b) /x = 0.5, for three values of 
molecule- lead coupling tL- U ' = 2, £m = 1, and T = 0.02. 



a quasi-particle peak at \i = 0. As discussed above, with 
decreasing tx,, a Coulomb pseudogap centered at \i = 
begins to be formed. As a consequence, at tz, = 0.5, 
ImEn(zo; n ) exhibits a finite value at small uo ni indicat- 
ing that electrons at the chemical potential have a finite 
lifetime inside the molecule at T = 0.02. Upon decreas- 
ing tjj further, ImEn(zo; n ) begins to approach the l/u n 
divergent behavior at small uj n , which corresponds to the 
limit of an isolated molecule. 

For comparison, panel (b) illustrates the self-energy at 
\i = 0.5, where the DOS exhibits a peak even at small 
values of t L (see -Gn(/3/2) in Fig. Efa) for t L < 0.5). 
In this case, ImEn(icj n ) remains approximately linear in 
uj n even at tL = 0.3, indicating that, outside the pseudo- 
gap region, the molecule maintains Fermi-liquid behav- 
ior. Eventually, of course, in the limit tL —> 0, the metal- 
lic behavior breaks down when the molecule no longer 
hybridizes with the leads. 



C. Kondo effect 

We now discuss the temperature dependence of the 
electronic structure of the N = 3 chain molecule. As 
shown in the Appendix, a cluster consisting of 5 lev- 
els adequately simulates a semi-infinite lead down to 
T « 1/800. Since the embedding potential is not af- 
fected by the size of the molecule, we use M = 5 lead 
clusters to investigate the molecular correlation effects in 
a wide range of temperatures. Fig. [7] shows the imagi- 
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FIG. 7: (Color online) Imaginary part of (a) interacting 
Green's function, ImGn (iuj n ), and (b) correlation-induced 
self-energy, ImEn(za; n ), for a linear molecule with N — 3 
for four inverse temperature values, /3 = 25, 50, 200, 500, and 
800. U = 2,t L = 0.5, t M = 1, and fi = 0.0 (half-filling). 



nary part of the {11} element of the interacting Green's 
function and the self-energy as a function of Matsub- 
ara frequency for T = 1/25, ... , 1/800 and tL = 0.5 at 
half-filling. The other parameters are the same as in 
Fig. [B^a). We note here that, in the limit of small uo ni 
— 7r _1 ImGn(icj n ) coincides with pi(uo = 0), the quasipar- 
ticle DOS at the chemical potential at site 1, as indicated 
in Eq. ([T4]) . As shown in panel (a), a sharp quasiparti- 
cle peak is formed at uj = when T decreases below 
about 1/500. Evidently, this peak may be identified as 
the Kondo resonance caused by the coupling between the 
localized spin in the half-filled second molecular level and 
the conduction electrons in the leads. The behavior of 
the self-energy is consistent with this trend, as shown in 
panel (b). While -ImEn(iu; n ) for T = 0.02 and 0.04 
extrapolates to a finite value in the limit of uo n — >• 0, at 
lower T it becomes linear in uj n . Thus, with decreasing 
T the system undergoes a transition from the Coulomb 
blockade regime to the Kondo regime. 

Since the second molecular level is energetically well 
separated from levels 1 and 3, its Kondo temperature 
may be estimated as follows. The wave function of this 
level in the site basis is given by 

|V2> = i(|l>-|3» 
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FIG. 8: (Color online) Partially integrated quasiparticle DOS, 
Pi(aO = —Gn(/3/2) (multiplied by /3), of a linear molecule 
with N — 3 as a function of chemical potential jjl for T — 1/50, 
1/200, 1/500, and 1/800. t L = 0.5, t M = 1, and U = 2. 

and the hybridization strength defined in Eq. (fl6j) is 

T 2 = 7Tti[p (0) + p 4 (0)]/2 = ^|po(0) = 0.131. 

Moreover, as pointed out above, the effective Coulomb 
interaction for this level is u = U22 = U/2 = 1. Thus, 
within the context of the half-filled Anderson model, we 
have u/A ^> 1, where A corresponds to T2 = 0.131. The 
Kondo temperature is then approximately given by the 
expression 

T K = (uA/2)^ 2 e- nu/8A ^ A/2u = 0.0157. (20) 

This estimate is fully consistent with the results in Fig. [71 
where both the Green's function and correlation-induced 
self-energy exhibit no noticeable changes as a function 
of temperature for T < Tk- Moreover, from the initial 
slope of Im En we obtain a quasiparticle weight Z « 1/8, 
yielding 53 Tk = ttZA/4 « 0.013, in reasonable agree- 
ment with the estimate quoted above. 

Figure [8] shows the temperature dependence of the par- 
tially integrated quasiparticle DOS at site 1, pi(fi) = 
— Gn(/3/2), for the same molecule as in Fig. [71 To com- 
pare different temperatures, the curves are multiplied by 
f3 since the weight function F(lj) in Eq. (fT5j) has inte- 
grated weight 7rT. As discussed above, the Kondo reso- 
nance appears at \i for T < T^, where Tk depends on 
fx. Apparently, Tk associated with the second molecu- 
lar level is larger than or comparable to the lowest value 
T = 1/800 in Fig. [8] Since the Kondo resonance appears 
at an energy close to fi = 0, it makes a large contribu- 
tion to pi(fi). As a consequence, the minimum (Coulomb 
pseudogap) between the two peaks for the second molec- 
ular level for T = 0.02 becomes shallower with decreasing 
T, and is eventually replaced by a broad single peak at 
T = 1/800. 

Interestingly, in contrast to this behavior of the second 
molecular level, the Coulomb pseudogaps for the first and 
third levels remain visible for the whole T range in Fig. [51 
The difference arises from the lower Tk values for these 



FIG. 9: (Color online) Partially integrated quasiparticle DOS, 
Pi(aO — — Gn(/3/2), of a linear molecule with N = 3 as a 
function of chemical potential \i for U = 0, 2, 4, and 6. £l = 1, 
t M 1, and T 0.02. 

levels. The wave function of the first level of the isolated 
molecule is 

|^i) = i(|l) + V / 2|2} + |3)), 

and the hybridization strength is 

ri = 7rt£[po(ei) + p 4 (ei)]/4 = 7Tt|p (ei)/2 = 0.046, 

where e\ denotes the centroid of the doublets, —\/2t — 
5C//16 = -2.04. By inserting A = T 1 and u = Uu = 
3/7/8 = 0.75 in Eq. Qgft for the case of half-filling, we 
obtain Tk = 2.4 x 10 -4 , which is much lower than the 
temperature range in Fig. [8l 



D. Large U region 

In Fig. O we fixed the onsite Coulomb repulsion U and 
varied the hopping integral between leads and molecule. 
Alternatively, it is of interest to inquire how the molecu- 
lar electronic structure changes when U is increased for 
fixed hopping. Here we consider the case of strong cou- 
pling where tL = tR = 1. Figure [9] shows the partially in- 
tegrated quasiparticle DOS at site 1, pi(fi) — — Gn(/3/2), 
for a chain molecule with N = 3 as a function of \i for four 
values of U. The other parameters are the same as in Fig. 
[51 As T/t L < 1, for U = 0, -Gn(P/2) is nearly identical 
to the non-interacting local DOS at site 1, pi(/i), except 
for a constant factor. The curve for U = 2 coincides with 
the one shown in Fig. [5fa). Compared with the bare 
non-interacting DOS, the three peaks are considerably 
broadened as a result of the intra- and inter-molecular 
orbital Coulomb terms appearing in Eq. (fT8|) . Moreover, 
the outer peaks are shifted to higher energies relative to 
the corresponding peaks at U = 0. When the Coulomb 
energy is increased to U = 4, all three DOS peaks be- 
gin to exhibit a minimum at their center. Finally, they 
evolve into double peak structures at U = 6. 
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Interestingly, the energy separations between the dou- 
ble peaks in Fig. [9] differ from those in the USCF ap- 
proximation discussed above (3/7/8 = 2.25 for the first 
and third orbitals and U/2 — 3 for the second orbital). 
This indicates that the off-diagonal Coulomb matrix el- 
ements ignored in Eq. (fT8j) become progressively more 
important with increasing U/t. Consequently, for [7 = 6, 
the three orbitals are significantly mixed by these off- 
diagonal terms. The spectrum may then more correctly 
be interpreted in terms of upper and lower Hubbard 
bands, each consisting of three peaks and split by the 
Mott-like gap at the center. Nevertheless, due to the 
proximity effect the DOS remains finite even at low en- 
ergies because the molecule is strongly coupled to the two 
metallic leads. 



E. Chain vs. ring molecules 

So far we have presented results for a chain molecule 
with N = 3. The results for other linear molecules that 
we have studied (TV = 2 to 5) are qualitatively similar and 
can be summarized as follows: (i) For a weakly correlated 
molecule, the quasiparticle DOS at the chemical poten- 
tial /i exhibits N peaks corresponding to the TV energy 
levels e m of the molecule, whose width is determined by 
the hopping integral between the lead and molecule, tz,. 
(ii) With increasing U/T m , these peaks start exhibiting 
minima at their centers, (iii) When U /T m is increased 
further, each quasiparticle DOS peak becomes a double 
peak structure, so that Pi(ja) = —Gu(/3/2) consists of 2N 
quasiparticle peaks as a function of /i, instead of the N 
peaks in the non-interacting limit. In the range of the 
correlation-induced pseudogaps, the electron self-energy 
exhibits a finite scattering rate. Furthermore, the elec- 
tron occupation of the molecule as a function of \i ex- 
hibits plateaus at odd integers, whose energy positions 
correspond to the location of the Coulomb pseudogaps. 
(iv) When temperatures is lowered to reach T^, which 
depends strongly on the Hamiltonian parameters and 
the molecular levels, Kondo resonances are formed in the 
Coulomb pseudogaps as a result of the strong coupling of 
the localized spin and conduction electrons in the leads. 

As an example, we plot in Fig. [TUla) the quantity 
—G av (P/2), which is defined as the average of Pi(fi) = 
—Gu(/3/2) over the TV molecular sites, for a linear 
molecule with N = 4 for two values of the molecule-lead 
coupling parameter, = 0.3 and 0.5. The other param- 
eters, U = 2, tM = 1, and T = 0.05, are the same as in 
Fig. [5j The comparison of Figs. 2] and [10] illustrates how 
the quasiparticle DOS peaks at the chemical potential 
\i evolve when the onsite Coulomb energy is increased 
from zero to a finite value. As mentioned above, there is 
a one-to-one correspondence between the DOS peaks in 
Fig. [Hand the double peak structures in Fig.[10l For the 
present molecule, the intra- and inter- molecular-orbital 
Coulomb energies are calculated as: U mrn = 3/7/10 
(m = 1 to 4), U 12 = U 13 = /7 24 = U34 = U/5, and 
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FIG. 10: (Color online) Partially integrated quasiparticle 
DOS averaged over all the sites, —G a v(/3/2), for a (a) lin- 
ear and (c) ring molecule with N = 4 for U — 2, £m = 1, and 
T — 0.02. Solid (blue) and dashed (green) lines correspond to 
th =0.5 and tz, = 0.3, respectively. The corresponding local 
densities of states for the non-interacting molecules (U = 0) 
are shown in Fig. [4[ Small numbers near the minima of these 
curves indicate chemical potentials giving integer occupations. 
Panels (b) and (d) provide the corresponding occupancies as 
a function of \i. 



U 1A = U 23 = 3/7/10. Thus, at positive /i, the DOS 
peaks of the third molecular orbital within the USCF 
approximation and in the small ti, limit are located at 
(VE - l)t/2 and (y/E - l)t/2 + 3/7/10, with a gap 3/7/10, 
while those of the fourth level appear at [y/E +l)t/2+C7/5 
and (\/E+l)t/2 + U/2, also with a gap 3/7/10. The peaks 
at negative \i are located symmetrically with respect to 
p = 0. The Coulomb pseudogaps in Fig. [TOTa) are in fair 



11 



agreement with these mean-field values. 

Finally, we discuss the ring molecule with N = 4 since 
it behaves quite differently from the corresponding linear 
molecule at low temperatures when it is half-filled. In 
Fig. [TUIb) we plot — G av {P/2) of this molecule for the 
same parameter set as for the linear molecule in Fig. 
[Tol a). With the molecular orbitals defined in Eq. (fT7|) . 
the intra-orbital Coulomb energies are U\\ — U44 = J7/4, 
U22 — U33 = U/2, while the inter-orbital ones are U mn = 
U/4 (m ^ n), except for U23 = 0. Hence, the DOS 
peaks of the first (fourth) level appear at —2t — U/2 and 
-2t - U/4 (2t + U/4 and 2t + U/2), with a gap U/4, 
while the degenerate peaks of the second and third levels 
are located at — U/4 and U/4, with a twice larger gap 
U/2. The peak energies in Fig. [TOl b) are in agreement 
with these USCF estimates. The ring molecule is half- 
filled when ji is located inside the pseudogap between 
the doublets at /i ~ ±C//4, with the first level essentially 
fully occupied and with the second and third ones singly 
occupied. 

The question arises as to whether the electrons in the 
second and third levels form a singlet or triplet state in 
the many-body ground stated"— which may arise as a 
consequence of off-diagonal Coulomb matrix elements ne- 
glected in Eq. JT5J). For the isolated molecule with the 
same U , we found that the singlet state has a lower en- 
ergy. Thus, because of the absence of a localized-spin 
degree of freedom, the N = 4 ring molecule at half-filling 
exhibits no Kondo proximity effect. To confirm this, we 
plot in Fig.[TT]the {11} element of the interacting Green's 
function and self-energy at half-filling. In striking con- 
trast to Fig. [7] for the linear molecule with N = 3, the 
correlation-induced self-energy is seen to preserve non- 
Fermi-liquid behavior at low temperatures, and p(uo = 0), 
i.e., the low frequency limit of — 7r _1 ImGii(zu; n ) in panel 
(a) does not exhibit a Kondo resonance. This explains 
why —G av (/3/2) in Fig. HDTb). when /i is located inside 
the pseudogap (\fi\ < C7/4), is much smaller than the cor- 
responding one for the linear molecule with N = 3 shown 
in Fig. [Ha), despite the fact that the pseudogap is nearly 
the same (~ U/2) for both molecules. 



IV. SUMMARY 

A new method for the evaluation of the electronic prop- 
erties of strongly correlated molecules coupled to semi- 
infinite metallic leads is proposed. By simulating the 
surface Green's functions of the leads in terms of small 
clusters, the many-body interactions of the combined sys- 
tem in the zero bias-voltage limit are obtained via ex- 
act diagonalization. The auxiliary energies and hopping 
terms of the lead clusters are derived by fitting the lead 
surface Green's functions at imaginary Matsubara fre- 
quencies. These fits are found to be sufficiently accurate 
to describe the Kondo physics, except in the limit of ex- 
tremely low temperatures. For moderate onsite Coulomb 
energies within the molecule, the density of states peaks 
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FIG. 11: (Color online) Imaginary part of (a) interacting 
Green's function, ImGn (i(jj n ), and (b) correlation-induced 
self-energy, ImEn (iuj n ), for a ring molecule with N = 4 for 
inverse temperature values, f3 — 50, 200, and 500. U — 2, 
th — 0.5, £m = 1, and \i — 0.0 (half- filling). 

of the non-interacting molecule are shown to split into 
doublets separated by correlation- induced pseudogaps. 
The molecular self-energy then exhibits a finite scatter- 
ing rate, as expected in the regime of Coulomb blockade. 
Outside the pseudogap regions, the self-energy retains 
ordinary Fermi-liquid behavior, characteristic of ballis- 
tic transport across the molecule. The one-electron hy- 
bridization between molecule and leads is shown to be a 
key parameter that governs the transition between the 
ballistic and Coulomb blockade regimes. If the chemi- 
cal potential is located inside a pseudogap, the molec- 
ular levels are integer occupied, so that a Kondo reso- 
nance appears upon lowering the temperature in the case 
of odd integer oocupancies. The present results suggest 
that the approach discussed in this work for molecules or 
quantum dots connected to metallic leads can describe, 
as a function of Hamiltonian parameters, electron fill- 
ing, and temperature, the full range of phenomena from 
Coulomb blockade to Kondo physics. In future appli- 
cations it would be interesting to apply this scheme to 
multiorbital dots and a variety of other models of inter- 
est for nanoscale devices. 
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FIG. 12: (Color online) Imaginary part of self-energy of 
single atom on semi-infinite lead along Matsubara axis at 
half-filling for temperatures corresponding to ft = 1/T = 
50, 100, 200, . . . , 1600 (from above), (a) U = 1 and (b) U = 2. 
£m = 1, th — 0.5. All results are obtained for M = 11, with 
a fixed Matsubara grid corresponding to /3m = 1600. 



Appendix: Single adatom on a semi-infinite lead 



FIG. 13: (Color online) Imaginary part of Green's function 
of single atom on semi- infinite lead. The parameters are the 
same as in Fig. 1121 The uncorrelated Green's function (U — 
0) is indicated by the green dashed curve. The symbols in 
panel (a) correspond to the Green's function l/(iuj n /Z + iA) 
for a single peak at /x = with Kondo temperature Tk = 
ttZA/4 (see text). 



In Section El we have discussed the main new fea- 
ture of the present scheme, namely, the simulation of the 
semi-infinite leads in terms of a finite set of levels. Es- 
sentially, the true embedding potentials which have con- 
tinuous spectra at real energies are replaced by those for 
finite clusters comprising a discrete set of poles. The cri- 
terion for this substitution is that along the Matsubara 
axis both versions of the embedding potentials should 
agree well for a given cluster size. Evidently, this fitting 
is accurate only at not too low temperature when the low- 
est Matsubara point is not too close to the real energy 
axis. Thus, for each cluster size, there should be a low 
temperature limit down to which the discrete set of levels 
accurately mimics the electronic properties of the actual 
semi-infinite lead. This limit may be determined by per- 
forming calculations for clusters with different sizes and 
by checking the consistency of the corresponding results. 

For this purpose we consider here the special case of 
a single atom (N = 1) adsorbed on a single semi-infinite 
lead which is equivalent to the single-impurity Ander- 
son model. The electronic structure of the correlated 
atom can therefore be directly compared with predictions 
within NRG. 1 — ^ If the lead is replaced by a cluster 
containing M levels, the calculation of the Green's func- 
tion of the combined system involves the eigenstates of a 
(M + l)-level cluster. To be specific, we choose hopping 
parameters t M = 1 and tL = 0.5. The Coulomb interac- 



tion is assumed to have the values U = 1 and U = 2. For 
the simple-cubic structure of the lead, the local density of 
states at \i = in the surface layer is p s (0) — 0.52/7T, so 
that the effective hybridization between atom and lead is 
A = 7rt L p s (0) = 0.13. In Kondo physics, it is customary 
to introduce the parameter x = U/(ttA) to character- 
ize the importance of Coulomb repulsion versus single- 
electron hopping. Thus, for U = 1 and U = 2, this 
parameter has the values x = 2.43 and x = 4.86, re- 
spectively. According to Eq. (|20]h the Kondo tempera- 
ture then has the values: Tk = 0.0157 for U = 1 and 
T K = 0.001 for U = 2. 

As discussed in Ref. |43|, finite-temperature exact diag- 
onalization can now be carried out for clusters involving 
up to about n s = 15 levels. Here we consider lead clus- 
ters up to M = 11, i.e., n s = 12. As shown below, 
these sizes are sufficient for temperatures down to about 
T = 1/1600, i.e., well within the Kondo regime for U = 1 
and above about 0.6 Tk for U = 2. Increasing M to 14 
would permit the study of even lower temperatures. 

Also, we point out here that, while the fitting of the 
lead Green's function in Eq. ([TO]) is usually done by us- 
ing the Matsubara points corresponding to the physical 
temperature T in Sec. IIII[ it is possible to introduce a 
fictitious Matsubara grid independently of T, which is 
used only for the purpose of fitting the lead Green's func- 
tion, i.e., for finding the parameter set {e^,^}. We de- 
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FIG. 14: (Color online) Upper (lower) panel: as in upper 
panel of Fig. [12] JT3j) except for M = 5. The Matsubara grid 
corresponds to /3m = 1600 for /? > 400 and to /3m = 400 for 
P < 400. The results for M = 11 at f3 = 1600 are indicated 
by the circles. 



note this fictitious Matsubara temperature by Tw Be- 
cause of the large cluster size (M = 11), we choose 
Tm = 1/1600, which should therefore provide excellent 
fits of goo in the entire range of physical temperatures 
considered, T = 1/1600, . . . , 1/50. The small value of T M 
implies that accurate low-energy behavior of the cluster 
Green's function Gf, is available. We caution, however, 
that, in using this technique, it is important to check that 
the cluster Green's function agrees with that of the lead 
not only on the fitting points in Eq. ([TP]) , but also on 
the real Matsubara points corresponding to the physical 
temperature. 

Figure [12] shows the low-energy behavior of the self- 
energy of the adatom for various temperatures. For 
U = 1 (upper panel) and T approximately less than 0.002 
(/3 > 500), Im£(zu; n ) is linear ~ iuj ni with a slope of 
about 6.4, yielding a quasi-particle weight Z « 0.135. At 
larger temperature, the self-energy develops a finite on- 



set, indicating a growing correlation induced low-energy 
scattering rate. The finite lifetime at /jl = then im- 
plies that the pinning condition of the interacting den- 
sity of states is increasingly violated. From the initial 
slope of the self-energy we can estimate the Kondo tem- 
perature by using the expression 53 Tk = 7rZA/4. Thus, 
Tk ~ 0.014 for U = 1, which agrees well with the esti- 
mate Tk = 0.0157 derived from Eq. ([20]) . According to 
the lower panel, for U = 2 the linear region of the self- 
energy is confined to much lower temperatures, with a 
quasi-particle weight of about Z < 0.01 and Tk ^ 0.001. 
This result is also in agreement with the estimate ob- 
tained from Eq. ([20]) . 

Figure [13] shows the quasi-particle Green's function for 
the correlated adatom in the same temperature range as 
in Fig.[H For U = 1 and T < 0.002, lmG(iuj n ) is seen to 
approach the same low-energy limit as the uncorrelated 
Green's function, as expected from the pinning condition. 
With increasing T, deviations from this condition become 
progressively larger, in correspondence with the behavior 
of the self-energy shown in Fig. [12] For U — 2, extrap- 
olation of ImG(i(jj — >• 0) to the pinning condition might 
be feasible only at the lowest temperature, in agreement 
with the estimate Tk ~ 0.001 given above. The behav- 
ior of G(iu) in the very-low-T region could be explored 
with greater accuracy by enlarging the lead cluster be- 
yond M = 11. 

To illustrate the accuracy of the adatom self-energy 
and Green's function for smaller lead clusters, we show in 
Fig.[Hthe results for U = 1 and M = 5. For f3 = 1600, a 
comparison with the results for M = 11 is also provided. 
Although the Green's function now is less accurate at 
low energies because of the reduced number of cluster 
levels, the self-energy still exhibits the correct qualita- 
tive trend: Below about T < 1/400, Im£(zco> n ) is lin- 
ear in cj n , whereas at larger T the low-energy scattering 
rate increases significantly. At low T, the initial slope 
of Im£(zco> n ) is nearly the same as for M = 11, yield- 
ing similar quasiparticle weight and Kondo temperature. 
Thus, in spite of the larger quantitative uncertainties in 
the case of the smaller lead cluster, the overall evolution, 
namely, from Fermi liquid behavior at low T to increas- 
ing low-energy scattering rates beyond about T = 1/400, 
is consistent with the more precise results for M = 11. 
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